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PREFACE 



The reasoning behind the simple, heat energy budget of a leaf is 
developed in this module into a detailed explanation of leaf transpira- 
tion. Although the mathematical model is conceptually simple, seven 
independent variables prevent easy exploration of the consequences of 
the model. The accompanying computer program provides the student with 
a tool to concentrate his efforts on the biology of heat and water balance 
of a leaf rather than time-consuming numerical calculations. Algebra 
and some knowledge of heat transfer physics are prerequisites. 
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INTRODUCTION 



Water is essential to life as we know it on earth. Water acts as 
a mediinn for carrying nutrients, ions, and suspended particles from the 
^oil into a plant and as a conducting system within a plant. Most organic 
^Issue is permeable to water and when the external water vapor pressure 
is less than that inside a plant, water will evaporate from it. Plant 
leaves have evolved stomates or pores through which to take in carbon 
dioxide from the surrounding air, but since the water molecule is of 
lighter molecular weight than the carbon dioxide molecule, it tends to 
escape readily from the leaf through the stomates. Normally the air in 
the intercellular spaces within a leaf is at or near saturation. This 
process of water evaporation from a plant is known as transpiration. 
Energy is required to convert liquid water to vapor. As a consequence 
of transpiration, energy is lost by a leaf and the temperature of the leaf 
is reduced. 

Temperature is extremely important to the physiological processes 
within a leaf including photosynthesis, respiration, growth, cell enlargement 
and division, protoplasmic streaming, and translocation. Most chemical 
and biological processes are temperature-dependent and as a general rule, 
they proceed slowly at relatively low temperatures, speed up in reaction 
with increased temperature, and at high temperatures become limited or 
terminated by adverse reactions such as the denaturation of proteins and 
the breakdown or dissociation of molecules. A typical biological response 
with temperature is Illustrated in Figure 1 for net photosynthesis. Many 
temperature zone plants have a maximum rate of photosynthesis at about 
25 to 30^C while plants of arctic or alpine habitats have an optimum 
temperature of about 15^C and some plants of extremely hot regions have 
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Figure 1. Net photosynthesle versus leaf temperature at high and 
low carbon dioxide concentration. From Lommen et al. 
1975. P. 41. 
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optima near 40 or 45**C (e.g., Lange et al. 1975, Bjorkman et al. 1972, 
Billings et al. 1971). When a plant of a cool habitat is functioning 
optimally at a particular temperature, a plant of a warm habitat will 
, scarcely function at all at that same temperature. Many plant leaves 
undergo irreversible thermal damage when their temperatures exceed 45 to 
48®C while some will not be damaged at 50 or even 55 ®C and others will be 
severely denatured at 42^C. Therefore, the temperature of a leaf is always 
of enormous significance to the viability of it. However, when one wishes 
to know the temperature of a plant or, in particular, the temperature of 
a leaf, it is important to ask why and with what accuracy. Is it necessary 
to know the leaf temperature to an accuracy of 0.1, 0.5, or 1**C, and why? 
The answer to this question will often determine the method of measurement 
and the precision necessary, or it will guide one with respect to the 
detail required in an energy budget calculation. 

The temperature of a leaf is the result, of energy flow between it 
and its surroundings. A leaf exchanges energy by the basic physical 
processes of radiation, conduction, convection, and evaporation, all of 
which go on simultaneously. If a leaf takes in more energy than it puts 
out» its temperature increases and if it loses more than it gains, the 
temperature will decrease. When energy in equals energy out, the leaf 
will ;naintain a constant temperature. 

Once the processes by which a leaf exchanges energy are recognized, 
it is straightforward to identify the properties of the environment which 
must be knoim* In fact, there are only two ways by which an environment 
interacts with an organism, including a plant leaf, and that is by energy 
flow or mass flow. In this module, we are dealing with both since the 
exchange of heat is an energy flow process and the diffusion of water vapor 
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is a gas exchange process which also involves the transport of energy. 

The remainder of the module is devoted to developing and analyzing 
simple models of how a leaf exchanges thermal energy with the environment. 
Although our knowledge of leaf morphology and physiology is more detailed 
than what is included in the models, the level of description here is 
sufficiently general to provide valuable insights. 

In the following pages, the reader is formally introduced to the 
leaf energy budget. Typical values for environment variables and leaf 
parameters are discussed. Some simple calculations are then made to see 
how radiation, convection and transpiration affect leaf temperature. 

The graphical method of analysis used by Gates (1968) is then 
employed to understand a more detailed energy budget model. Problems 
with accompanying solutions are followed by an appendix on symbols, units 
and dimensions. Finally, a second appendix outlines how the student could 
use a computer program called TRANS to answer problems and explore conse- 
quences of the second model. 

LEAF ENERGY BUDGET 

Radiation incident upon a leaf includes shortwave radiation in the 
form of direct sunlight, skylight and reflected light, and longwave 
radiation emitted by the atmosphere and clouds, or by the surfaces of 
ground, nearby plants, and other objects. 

A plant may absorb AO to 80 percent of incident sunlight and skylight 
and it will absorb 96 percent of incident longwave thermal radiation. 
A plant leaf redistributes this absorbed radiation by emitting a substantial 
amount of it as longwave radiation from the leaf surface. The remainder 
is partitioned into convective and evaporative heat exchange and a very 
small amount is utilized for photosynthesis. Only certain wavelengths of 
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light are effective for photosynthesis. The amount of energy converted 
to biomass is very small (a few percents of the incident sunlight at most) 
and may be neglected in considerations of the energy budget of a leaf. 
The amount of radiation emitted by a leaf is proportional to the fourth 
power of the absolute temperature of the leaf surface. This phenomeno- 
logical relation is known as the black body radiation law and is written 
R = a[T+273]\ 

The flow of heat by convection is proportional to the temperature 
difference between the leaf and the air. Convective heat transfer is 
proportional to the wind speed and varies inversely with the characteristic 
dimension of the leaf which affects the boundary layer thickness. Because 
of the viscosity of air, there is an air layer which adheres to any surface, 
including that of a leaf known as a "boundary layer." When the wind blows, 
the boundary layer of viscous air represents a transition zone between 
zero air flow at the leaf surface to free air flow at some distance from 
the leaf. A temperature difference may exist across a boundary layer of 
air and it is across this that heat moves from the leaf to the air by 
conduction and convection. The larger a leaf, the thicker the boundary 
layer of air adhering to the surface. A small leaf will have a boundary 
layer about 1 mm thick and a large leaf, like a banana leaf, will have 
a boundary layer of 1 or 2 cm thickness. It happens that for broad flat 
leaves the rate of heat transfer by convection is proportional to the 
square root of the ratio of wind speed to leaf dimension. Increase wind 
speed fourfold and the rate of convective heat exchange will increase 
twofold or increase leaf dimension fourfold and convection will decrease 
twofold. The characteristic dimension is approximately the average width 
of the leaf. Although the length of a leaf affects convection in the 
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interest of simplicity, a single characteristic dimension given by the leaf 
width is used here. The more complex relationship involving leaf width and 
length is introduced later. 

The transpiration rate is determined by the water vapor pressure 
or density difference between leaf and air. This relationship is described 
in more detail in the next two sections. Energy is required to convert 
liquid water to vapor and the amount of energy is known as the latent heat 
of water. It is temperature-dependent but at 30^C its value is 2. A3 x 10^ 
J kg""^. Each of the environmental variables affect the leaf energy status 
simultaneously. They are radiation, air temperature, wind speed, and water 
vapor pressure, density, or relative humidity. A plant leaf responds to 
these by assuming a certain temperature and by losing water at a particular 
rate depending upon the properties of the leaf. Once the energy budget 
relationship is written, the student will see how it is that each of these 
environmental variables enter the common energy pool for the leaf. 

The steady state energy budget for a leaf is 

Thermal energy gained = Thermal energy lost 

Absorbed radiation = Reradiation + Convection + Evaporation 

h 



Q = ea T + 273 ^ + k 
a x> 



[T^ - T^] + L(T^)E (1) 



where 



0 ■= total amount of radiation absorbed in W m , 
^a 

e e emissivity of leaf surface to longwave radiation, 

a « (Stefan-Boltzmann constant) 5.67 x 10""® W m""^ K~**j 

T^ « leaf temperature in ^C, 

kj « 9.14 J m-^s-^ ^C-^ 

V ».wind speed in m s~^. 
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D » leaf width in m, 

T " air temperature in **C, 
a 

L(Tj^) » latent heat of vaporization of water in J kg"^ as a function 
of leaf temperature and is equal to 2. A3 x 10^ at 30''C and 
2.50 X 10* at O^C, 

E « transpiration rate in kg m"^ x s"^ , 

For any given value of E, a specific value of T^^ will balance Equation (1), 
provided all other quantities are known. However, E depends on a unique 
set of environmental variables. It is seen in the next paragraphs that 
E is a function of T.. Therefore, it is possible to determine E and T« 
simultaneously. If the air temperature is warmer than the leaf temperature, 
heat is gained by convection rather than lost and this term becomes negative 
on the right-hand side of the energy budget equation. 

Resistance to Water Loss 
There are three requirements for the loss of water from a leaf. 
There must be water available in the leaf, there must be energy available 
to convert liquid water to vapor, and finally, there must be a vapor 
pressure of density gradient along which water vapor may flow from inside 
to outside the leaf beyond the bounday layer of air which adheres to the 
leaf surface. Liquid water at the mesophyll cell walls within the leaf 
is vaporized and from the intercellular spaces, this water vapor passes 
out of the leaf by diffusion through the stomates or through the leaf 
cuticle. Usually, the cuticle is coated with a wax layer which is 
relatively impervious to water and as a result most of the water lost from 
a leaf is through the stomates. As with fluid passing through any tube or 
pipe, there is resistance to vapor flow by viscous drag with the walls. 
The resistance to water vapor diffusion through the substomatal and 
Btomatal spaces with the leaf is rp in s m"^ and represents an average 
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value for the entire leaf. In addition, the water vapor must diffuse 

across a boundary layer of air adhering to the leaf surface and this offers 

a resistance r^ given in s m"^. The boundary layer resistance is affected 

^ by the wind speed across the leaf. In fact, the greater the wind speed 

^the thinner the boundary layer thickness; and hence, diffusion resistance 

is greater on large leaves than on small leaves. Experimental results show 

that the boundary layer resistance varies directly as the square root of 

the leaf dimension (width) D and inversely with the square root of the wind 

speed V, a relationship which is precisely the inverse of the influence of 

wind speed and leaf dimension on the exchange of heat by convection. 

Therefore, the boundary layer resistance r is given by 

a 



^^8 = ^2 



(2) 



Although there are complexities concerning air flow about a leaf which 
cause kg to change its value for leaves of characteristic dimension less 
than 0.005 m, we shall use a single value here for all leaf dimensions, 
kg « 200 s^'^^m'*^ when D is in m and V is m s"^. Since the ratio of D/V is 
involved here, nothing is changed if D is in cm and V in cm s""^. 

Transpiration Rate 

The rate at which water vapor will diffuse out of a leaf will 
depend directly on the difference in water vapor density inside and outside 
the leaf, just as the rate at which water will flow downhill depends on the 
difference in height between the top and bottom of the hill. 

For the moment, it is assumed that the air in the intercellular air 
spaces inside the leaf is at saturation at the temperature of the leaf and 
has a saturation density gd^(T^). The air outside the leaf and beyond. the 

boundary layer has a water vapor density eq\t^l>to [r.h.] x d (T ), where 

1 i'*f s a a 
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r.h. is the relative humidity and d (T ) is the saturation water vapor 
density of the air as a function of the air temperature. The relative 
humidity is, by definition, the ratio of the actual water vapor' density 
of the air to the saturation density. Usually, it is given in percent, 
But as used here, it is a 'ecimal fraction. The rate of water loss from a 
leaf per unit area per unit time is equal to the gradient for water vapor 
density divided by the resistance to water vapor movement and is given by 

d (T ) - [r.h.] d (T ) 
a ^a 

where r is given by Equation (2). Here, if d (Tp) and d (T ) are given 

** s X/ X/ s a a 

in kg m"^ and r^^ and r^ in s m"^, then E is in kg m"^ s"^. 

COMPLETE ENERGY BUDGET 
By substituting Equation (3) into Equation (1), and including the 
expression given by Equation (2) in Equation (3), one gets for the full 
expression of the leaf energy budget the following: 

h 



% " ^^^'^l 273]** + kj 



V 



- ^a^ 



d (T ) - [r.h.] d (T ) 

+ L(T^) ^ ^ (4) 



Here all environmental variables which affect the energy status of the leaf 

act simultaneously. They are Q , T , e, D, and r., while the absorptivity 

a a jC 

of the leaf surface to radiation is buried in the term Q . When all the 

a 

environmental variables and appropriate leaf properties are known, a unique 
value of Tj^ will balance this equation. 
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Values of the Leaf Parameters 

Leaves are of many sizes from 1,0 x 10"^ m by 2.0 x 10""^ m for a 
Douglas fir needle to 0.3 m by 1.5 m for a banana leaf. Some leaves will 
have lengths very nearly equal to their widths. Internal diffusion resis- 
tances of leaves vary from less than 100 s m""^ to infinity but with commonly 
occurring values between 200 and 2000 s m""^. Most leaves will have an 
absorptivity to shortwave radiation of about 0.6, but it may vary from 
0.4 to 0.8. Longwave absorptivity will equal the longwave emissivity of 
0.96. The total amount of radiation absorbed by a leaf will vary from about 
400 to 800 W m""^ but higher and lower amounts may be encountered. 

Values of the Environmental Variables 

It is relatively easy for any of us to visualize the range of values 
for air temperature. A leaf is likely to encounter freezing when the air 
temperature is 0°C. Maximum air temperature for hot desert conditions, for 
example, would be about 45^C. Relative hxjmidity is also easy to visualize 
as a result of common experience. For example, humid air would have a 
relative humidity of about 70% or greater and for very dry air r.h. «= 30% 
or less. It is easy for one to look up in tables the water vapor density 
of saturated air at any particular temperature. Values range from 1.289 
kg m""^ at O'^C to 1.068 kg m"^ at 45°C. Wind speeds are a part of common 
experience. Still air usually has very slight air movement which we 
arbitrarily put at 0.1 m s"^. A gentle breeze is 1.0 m s"^ (2.2 mph) and 
a moderate wind is 10 m s~^ (22 mph). Radiation is the least familiar of 
all of the environmental variables and yet it is the most significant of them. 

INFLUENCE OF ENERGY COMPONENTS ON LEAF TEMPERATURE 
It is useful to evaluate the major components separately of the 
leaf energy budget in order to get a feeling for their values and influence. 
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If a non-transpiring leaf is ina vacuum, its temperature is determined 
solely by radiative exchange and its energy budget is 

- eo[Tj^ + 273] \ (5) 

,Let e - 0.96. If - 800 W wT^ , T^^ - 75.2'*C; 600 W m'^ , T^^ - SLl-'C; 
and - 400 W m~^, T^^^ " 19.8''c. A leaf tein)erature of 75''C is very warm 
indeed and no plants as we know them would survive such a temperature. 
Even at 51 most plants would sustain serious heat damage. However, 
plants live in air, a, fluid, which flows about the plant and takes away 
excess heat by convection. Just as a plant loses heat by radiation, it 
must lose heat by convection as well, and this it cannot avoid. We shall 
now determine the extent to which a plant leaf may have its temperature 
influenced by convection. 

If the non-transpiring leaf is placed in air, the leaf temperature 
is significantly affected by the air temperature. If the ratio V/D is 
large, convection becomes very significant and the leaf temperature is 
tightly coupled to the air temperature.^ Large V/D suggests high wind 
speeds or small leaf dimension. If V/D is small, the leaf temperature 
is less strongly affected by the air temperature, but still it is signifi-- 
cantly influenced. The energy budget is given by 
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If the air temperature is 30*^0, one can determine the exact 
numerical relationship between and Tj^ for any given ratio V/D. For 

example, we will determine values of T^. for Q « 800, 600, and 400 W m 

x» a 

for V/D « 1.0, 10.0, and 100. These values are listed in Table I. The 

easiest way to solve Equation (6) is to put in values of V/D, T^,, and T , 

x» a 
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Radiation 
. only 




Radiation and 
convection 
V/D(8-^) 



Radiation, convection, 
and transpiration 
V/D(8-^) . 



1.0 


10.0 


100.0 


1.0 


10.0 


100.0 


61.5 


39.7 


33.5 


36.8 


31.8 


30.3 


39.1 


34.0 


31.2 


30.4 


28.5 


28.7 


26.1 


28.3 


29.4 


23.0 


24.9 


27.1 



Table 1. Leaf temperature In for various amounts cf absorbed 
radiation exchanged by radiation only, by radiation and 
convection, and by radiation, convection, and transpiration. 
Values of variables used are Ta«30^C, r.h.«50%, and r^ « 
100 s m""^. The addition of convectlve exchange to the 
energy budget will decrease the difference between leaf 
temperature and air temperature compared to the leaf 
temperature computed when only radiation is included. 
As the ratio V/D increases, the magnitude of the difference 
between leaf temperature and air temperature decreases. 
Note that if absorbed radiation levels are above black body 
value (for Qg-SOO and 600 W m"^ because here black body level 
is Qa»CTT^« 5. 67x10-® x (273+30)^ - 478 W m"^), T^-Tg is 
positive and if Q3 is below the black body value T^-Tg is 
negative. If transpiration is included, it will reduce all 
leaf temperatures from the values computed when only 
radiation and convection transfer are included* 



EKLC 



17 



-13- 

determlne and then plot Tj^^ versus at fixed ratios of V/D and a fixed 

value of T ; in this case T » 30^C. 
a a 

For V/D « 1.0, simply placing the non- transpiring leaf In air of 

temperature 30^C dropped the leaf temperature from 75.2 to 6l.5**C at 

■ 800 W m"" , a decrease of 13.7^C. Clearly, convection may play a very 

strong role in influencing leaf temperature. If V/D « 10.0, the decrease 

is, of course, even more dramatic. At Q ■ 800 W m"^ and T » 30**C, 1^ » 

a a Jc 

39.7°C, there is a decrease of SS-S^C from the pure radiation regime. 
If V/D = 100 and Q -= 800 W m"^, the leaf temperature is 33.5°C with 
convection, a change of 41.7®C. 

We tend to forget the significance of convection in the ordinary 
world but when we make comparisons of this kind to show what would happen 
to leaf temperature without convective cooling, the influence of convection 
becomes very clear. At lower amounts of absorbed radiation, the influence 
of convection is less strong, but still important until the leaf temperature 
is equal to the air temperature. Then, when the quantity of absorbed 
radiation is very low, which will occur for an exposed leaf on a clear 
night, the leaf temperature is below the air temperature and convection 
will deliver heat to the leaf, thereby, warming it. This is seen in Table I 
for = 400 W m~^; when, without convection, the non-transpiring leaf had 
a temperature of 19. 8*^0, and the in air at 30® C, the leaf temperature 
increased to 26.1, 28.3, and 29.4''C, respectively, at V/D « 1.0, 10.0, 
and 100.0. 

The influence of transpiration is now readily seen by including the 
transpiration term in the full energy budget expression as given in 
Equation (4). Now one must decide on a value for the internal resistance 
of a leaf and the relative humidity of the air. Let r^ « 100 s m""^ and 
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r.h. 50%. It is more difficult to solve Equation (4) then it is to solve 
Equations (5) or (6). Graphical methods are possible although cumbersome. 
It is far easier to use a calculator or a computer and this is now done. 
The values are listed in Table I. 

At Q » 800 W m-S T = 36.8, 31.8, or 30.3''C for V/D « 1.0, 10.0, 

^ a x> 

and 100.0, respectively. Hence, transpiration has reduced the leaf 
temperature by 24.7, 7.9, and 3.2^C for the three ratios of V/D, as shown 
by comparing the values listed in columns 6, 7, and 8 with those in columns 
3, 4, and 5. These temperature reductions caused by transpiration are quite 
striking and would be significant in terms of photosynthetic and respiration 
rates. At lower levels of absorbed radiation, the influence of transpiration 
if reduced, since there is less energy available for the evaporation of 
water. Transpiration always produces leaf temperatures lower than the case 
of radiation and convection only. 

MORE DETAILED ENERGY BUDGET 
So far, the energy budget equation has contained only the leaf 
width as the characteristic dimension for convective heat transfer and 
gas diffusion. Such a simple expression is only strictly correct for leaves 
which are nearly square; however, realtively few leaves are of this form. 
Most leaves have a long dimension which is considerably greater than the 
narrow dimension. Some years ago, Gates (1968) reported investigations 
of the heat and water vapor transfer from simulated leaves of blotting 
paper of a variety of lengths and widths. He found if D is the characteristic 
dimension of the leaf in the direction of air flow and W is the dimension 
transverse to the wind, that the energy budget is given by 

o 
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% - ea[T^ + 273]^ + i^f [T^ - Tj 




,0.50 



(7) 



- where 



kj « 9.14 J m- 



,-2 g_}s o^-i 



and 



- 183 s"^ m-^ 



Gates elected to include W in the evaporation term only since the 



amount of data was not sufficient to reveal the functional relationship 
in both the convective and evaporative terms simultaneously. Equation (7) 
is more realistic than Equation (A), which, nevertheless, is a good approxi- 
mation. If the wind is blowing across the leaf in the direction of the 
leaf width, then, the characteristic dimension D is the average leaf width 
and W is the average length. If, on the other hand, the wind is blowing 
along the length of the leaf, then, D is the average length and W is the 
average leaf width. 



variables, leaf temperatures and transpiration rate, as a function of any 
conceivable combination of values of the independent variables and leaf 
parameters which appear in Equation (7). The four independent variables 
used are radiation absorbed, air temperature, relative humidity of the air, 
and wind speed; and the three leaf parameters are characteristic dimension, 
leaf width and internal resistance. Gates and Papian (1971) published a 
tabulation of leaf temperatures and transpiration rates for various sets 
of values of the dependent variables and leaf parameters. Similar plots are 
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Computer Leaf Temperatures and Transpiration 



It is straightforward to obtain values of the two dependent 
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shown in a paper by Gates (1968). These plots have proved useful in under- 
standing Equation (7), and therefore at this time we will undertake an 
explanation of how they are produced. 

Let us begin by specifying Q « 800 W m"^, T » AO^C, r,h. ■ 0.20, 
- a a 

D « 0.05 m, W «= 0.05 m, r «= 0.0, s m"^ and V >= 0,10 m s"^. Now all the 
5 a 

parameter values of Equation (7) are known except leaf temperature. Using 
numerical techniques or a trial-and-error procedure one can find T^. Once 

is known, the transpiration rate is calculated using the last term of 
Equation (7). It is then possible to plot this point (as a triangle) in 
Figure 2a. If V is changed to 2.1 m s"^ and the other variables are held 
constant, we get another pair of transpiration and leaf temperature points 
to plot shown by a circle in Figure 2a. If V = A.l m s"^ and 6.1 m s"\ 
the square and star points respectively can be computed in the same 
fashion. 

These four points are then connected as in Figure 2a with l^s. 
This line describes all the pairs of leaf temperature and transpiration rate 
points with the fixed conditions specified above and wind speed varying from 
0.10 m s"^ to 6.1 m s"\ 

It is also possible to fix wind speed and let leaf resistance vary. 
As an example, assume the same conditions as at the beginning of the previous 
paragraph but let V = 2.1 m s"^ and r^ vary continuously from 0.0 to 
600 s m"^. The result is the line of 6's shown in Figure 2b. Using this 
procedure, a series of eight lines is generated for each plot. If the plots 
look confusing, simply locate a point of intersection of two lines and ask 
what are the parameter values. It is then easy to find out which parameter 
is varying for any given line. 
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Figure 2a. Plot of transpiration versus air 
temperature using the energy balance equation 
(7) from the text, with Qj,=800 w m^. Tg=40''C, 
D=0.05 m, W=0.05 m, and rfO.O s r\ For 
V=0.10, 2.1, 4.1 and 6.1, m s"' symbols A, 
°, p and * are used. The line of I's is 
of V varying continuously from 0.10 to 6.1 
m s"'. 



Figure 2b. Same as 2a except that the line 
of 6's is for V=2.1 m s~^ and r^ varying from 
0.0 to 600 m s"^ 
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Sample Plots of Transpiration and Leaf Temperature 
It Is possible to select as changing variables in the graph any two 
of the following Independent variables and leaf parameters: amount of 
.^radiation absorbed, air temperature, wind speed, relative humidity, leaf 
^dimensions, and diffusion resistance. For our first example, we shall use, 
as changing variables, wind speed and internal diffusion resistance. Let 
wind speed values by 0.1, 2.1, 4,1, and 6.1 m s~^ and Internal diffusion 
resistance values be 0, 200, 400, and 600 s m-^ 

A leaf size of 0.05 x 0.05 m is selected for the first examples. 
Plots are generated for the following conditions: 
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r .h. 


- 0.20 
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T " 
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lO^C; 


r.h. 


- 0.50. 



Plot 3a has many interesting features. It is characteristic of a 
hot, dry, high intensity radiation environment, such as one finds in a 
desert. In still air, represented by V « 0.1 m s"^, a leaf with a high 
diffusion resistance becomes very warm. For r^. 600 s m"^, T^ « 45^C, 
but if = «, Tj^ is nearly 15''C above the air temperature of 40''C. 
This can only be obtained by extrapolating the straight line for V « 0.1 
m s""^ to where it intersects the temperature axis at zero transpiration. 
(The 5's line. It would be easiest to see this if one enlarged the scale 
of plot and let the ordinate begin at 0.0 transpiration rate.) It is 
evident that all the constant wind speed lines intersect and cross over at 
Tj^ « T^ « 40*^0. Note the small decrease in transpiration rate for 

■ 600 s m~^ when T^^ < 42 ®C and the wind speed has increased to V = 2.1, 
4,1, and 6.1 m s~^. This is an interesting phenomenon and should be 
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Plot 3a. Leaf transpiration versus leaf 
temperature with leaf resistance rj^ and wind 
velocity V varying. 



Plot 3b, Enlarged 3a. 
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understood by the reader. Lopk at Eq[uation (7) and consider the effect of 
increased wind speed. An increase of V, when the leaf temperature is above 
the air temperature, produces a decrease of leaf temperature by convective 
cooling. The increased wind speed also diminishes the thickness and 
^resistance of the leaf boundary layer. This, by itself, would cause the 
'transpiration rate to increase. However, the convective cooling of the 
leaf caused the water vapor pressure inside the leaf to drop and that 
produces a reduction of transpiration. In effect, what happens is that 
there is a greater reduction of the numerator in the evaporation term of 
Equation (7) (the driving force) than there is of the denominator (the 
resistive force). This phenomenon shows up at relatively high values of 
internal leaf resistance. At lower leaf resistance values, of AOO s m""^, 
when the wind speed increases the transpiration rate increases and the line 
of constant internal leaf resistance turns upward with wind speeds greater 
than about 1 m s"^ (see Figure 3b, line 3)'. 

Now consider what happens when the leaf diffusion resistance is 
small, e.g., 200 sm"^^. Here in Plot 3b, the line at constant r^ « 200 s m-^ 
is doubled valued in E for certain T^^^. The line curves back on itself as 
wind speed increases. At firsts as the wind speed increases from 0.1 m s""^ 
to about 1.0 m s"\ the leaf temperature decreases as the transpiration rate 
Increases, but at higher wind speeds; although the transpiration rate 
continues to increase, the leaf temperature begins to increase. At very 
high wind speeds, it would approach the air temperature because convection 
is transferring heat to the leaf Bince the air is warmer than the leaf. 
A general feature of Plot 3b is that for hot, dry, high radiation conditions, 
a leaf is warmer than the air temperature if the internal diffusion resis- 
tance is above a certain value-^here, estimated at about 300 s m""^* If the 
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leaf resistance is less than this value, the leaf temperature is below the 
air temperature. 

Plot 4 represents leaf temperature and transpiration rates for 
.warm, dry, moderate radiation conditions, with V and r^ as changing 
cjuantities. Only leaves with relatively low internal resistance have 
temperature less than the air temperature (line of I's). Leaves with 
resistances of about 200 s m*^ or greater will be from 0 to about S'^C above 
air temperature, depending on the wind speed (lines of 2's, 3's, and 4's, 
Plot 4b). Where does the line of zero internal diffusion resistance show 
up in this plot or is it beyond the scale used? 

Plot 5 is for warm, humid, moderate radiation intensity conditions 
such as occur in many tropical situations and during moist, summer temperate 
region days. Transpiration rates are generally low, since the vapor 
density difference between leaf and air is not veiy great. Now the leaf 
temperatures are always above the air temperature for all realistic values 
of r^. 

Plot 6 represents cool, moderately humid, low radiation conditions 
and exhibits some interesting features. All leaf temperatures are below 
the air temperature. An increase of wind speed always warms the leaf and, 
at constant internal leaf resistance, produces a small increase of trans- 
piration. Generally, the transpiration rates are very low because of the 
limited amounts of energy available. The line of zero leaf resistance 
constant is well displayed here. The reader should ask if the conditions 
used here, e«g., Q « 300 W m-^ T^ « lO'^C, and r.h, « 0.50 are realistic 
and, if so, when and where will they occur. 

Now it is possible to explore other dimensions of the transpiration- 
leaf temperature plot. Let the changing variables be air temperature and 
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Plot 4a. Leaf transpiration versus leaf 
temperature with leaf resistance rjj and wind 
velocity V varying. 



Plot 4b. 4a enlarged. 
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Plot 5, Leaf transpiration versus leaf 
temperature with leaf resistance rj^ and wind 
velocity V varying, 



Plot 6. Leaf transpiration versus leaf 
temperature with leaf resistance rj^ and wind 
velocity V varying. 
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internal diffusion resistance. Use air temperature of 10, 20, 30, and 

AO^'C and resistance of 0. 333, 667, and 1000 s m'"^. The leaf size is 

0.05 X 0.05 m as before. Plots are generated for the following' conditions : 



7) 




700 W 




V -= 


0.1 


m 




r.h. 




0.20 


8) 




700 W 




V » 


0.1 


m 


B-'i 


r.h. 




0.80 


9) 


Qa = 


700 W 




V = 


2.0 


m 




r.h. 




0.50 


10) 


^a" 


300 W 


in-2; 


V " 


0.1 


m 




r.h. 




0.50. 



Plot 7 represents an environment of no wind, moderate radiation 

intensity, and very dry air. There are lines of constant air temperature. 

Note where it is in the plot that T. > T and T. < T . Locate the transition 

a a 

point on the air temperature lines and draw a light line connecting these. 

Plot 8 is for an environment of no wind, moderate radiation intensity 

and high humidity. Compare this with Plot 7. Notice the change in slope of 

the constant temperature lines and the curvature of the constant resistance 

lines. It is seen there is an optimum temperature for the freely evaporating 

surface with zero internal resistance. Why? As the air temperature increase 

from 10 to. 20 to 30^C along a constant resistance line, the leaf temperature 

is increasing but at the same time the difference T^^^ - T^ is getting smaller. 

Since the term for radiation emitted by the leaf increases with the fourth 

power of T^, it is taking up a disproportionate amount of the available 

energy which is constant at Q « 700 W m~^. Not only is there less energy 

a 

available for convection, but less for evaporation of water. Also, as T 

a 

increases, the vapor density difference in the evaporation term of Equation 
(7) becomes less and the driving force to remove moisture from the leaf 
diminishes. 

Plot 9 is for conditions of moderate radiation intensity, moderate 
humidity » and light wind. The rate of transpiration has increased over that 
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Plot 7. Leaf transpiration versus leaf 
temperature with leaf resistance r^^ and air 
temperature Tg varying. 



Plot 8. Leaf transpiration versus leaf 
temperature with leaf resistance and air 



temperature T varying. 
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ot 9. Leaf transpiration versus leaf 
mperature with leaf resistance r^ and air 
mperature T varying. 



Plot 10. Leaf transpiration versus leaf 
temperature with leaf resistance r^ and air 
temperature T varying. 
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in Plots 7 and 8. Here the leaf temperature is above air temperature for 
high resistance values, but for low resistances it is below; whereas in 
Plot 8, it was above air temperature for nearly all of the graph. 

Plot 10 is for a low amount of absorbed radiation, no wind, and 
moderate humidity. The constant air temperature lines are very steep and 
the leaf temperature is always less than the air temperature. Note that 
for each constant resistance line, there is an optimum air temperature for 
which transpiration is a maximum. Why does transpiration diminish at high 
air temperatures? Are the environmental conditions realistic throughout 
the graph? Is the amount of radiation absorbed by the leaf consistent with 
the air temperature regime throughout the plot? If not, draw a line on the 
graph demarking a region of real environments from unreal ones. 



In this module, two simple models of the thermal energy budget of a 



leaf have been introduced. Radiation, convection and evaporation, the most 
important heat transfer processes, are included in the models (note that 
conduction of heat along the leaf stem is not included) . Because of the 
multidimensionality of the models, a graphical analysis was used. The 
models show that 



1. leaf temperature can be above or below air temperature, 

2. the larger the difference between absorbed radiation and the 
black body level of radiation, the larger the difference between 
leaf temperature and air temperature, 

3. the greater the wind speed, the smaller is the difference 
between leaf temperature and wind speed, and 

4. Increasing evaporative water loss via transpiration always 
decreases leaf temperature « 

In this formulation, the plant "can control" leaf temperature or 



transpiration by its size (D and W) and by its physiology (k, and r.). 



CONCLUSION 
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There are, however, a nimber of leaf characteristics which have not 
been explicitly or Implicitly included in the leaf energy budget. These 
include stomata size and shape, wax cuticle properties, leaf orientation 
and leaf absorptivity to shortwave radiation^ 

Finally this work has help lead to research in related areas. 
Scientists have developed equations to model (1) photosynthesis which couples 
the interchange of water vapor and CO^ (Lommen et al. 1975; Tenhunen, Yocum 
and Gates 1976; Tenhunen et al. 1976), (2) whole-plant water transport 
(Famum 1977), and (3) optimal leaf form (Taylor 1975; Parkhurst and Loucks 
1972; Givnish and Vermeij 1976), 
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PROBLEMS 

1. Using the example of a leaf vith dimensions .05 m by .05 m, Dr. Gates 
shows in his paper that under certain circumstances an increase in 
wind speed will increase transpiration while in other cases an increase 
in wind speed will decrease transpiration. Find examples of these two 
cases for a Douglas fir needle with D » .001 m, W « .02 m. 

2. Your light plane has just been forced down at 5000 feet in the Andes 
mountains. It is a hot (SS^'C), still (wind speed - 0.1 m sec"^), and 
dry (r.h. » 20 percent) day in late summer. To pass the time until 

you are rescued, you look at the surrounding flora and become interested 
in leaf sizes. The sun is shining brightly and you guess that, without 
as much atmosphere to absorb the radiation, the plant leaves must be 
absorbing about 1045.6 W m""^ of radiant energ , You also note that the 
ground is dry so you surmise that the plants are water-stressed and must 
have internal diffusive resistances of about 2000 s m"^. Assuming that 
these are the conditions which determine plant survival in this area, 
and assuming that plant leaves cannot tolerate temperatures in excess 
of A5^C for extended periods, what are the shape and size of the 
largest leaves you expect to find? 

3. Equation (7) in the module indicates how a leaf dissipates the radiant 
energy it absorbs. The three terms on the right-hand side of (7) 
correspond to reradiation; loss of sensible heat by convection; and loss 
of latent heat through transpiration. For each of the three situations 
listed below, find the relative contribution of each of these mechanisms 
to the dissipation of the Incident radiant ener';:y. 

(a) In Figure 3, the situation indicate^ by u . • intersection of the 
RL " 200 sec m""^ and V « 0.1 m .^e lines, 



40 



-30- 



(b) In Figure 3, the situation indicated by the intersection of the 
RL ■ 0 sec m"^ and B » 2 m sec"^ lines. 

(c) In Figure 8, the situation indicated by the intersection of the 
RL » 5000 sec m"^ and TA - O^^C lines. 
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PROBLEM SOLUTIONS 

1. Transpiration Increases with wind speed for lines 1 and 2, J..e., 
RL - 168 s n-*, but decreases with wind speed for lines 3 or 4 
Csee Plots 11a and lib). 

2. The leaves will be long and thin so that their dimension In the 
direction of the wind Is always small. They must be less than about 
.007 m wide, and are essentially unlimited in length (see Plot 12). 

3. 
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Plot 11a. Leaf transpiration versus leaf 
temperature for Problem 1. 



Plot lib. Leaf transpiration versus leaf 
temperature. The lover right-hand corner of 
Plot A is blown up to show lines 2, 3, and 4. 
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Plot 12. Leaf transpiration versus leaf temperature for Problem 2. 
Diameters above 0.0067 m make the leaf temperature raise 
above 46''C. 
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TRANSPIRATION AND LEAF 


TEMPERATURE 






APPENDIX 1. Symbols, Units 


, and Dimensions 




Symbol 


Quantity 


Units 


D imen s i on 


D 


Leaf dimension parallel to the 
wind direction 


M 


L 


d 
s a 


Saturation density of the air 


kg m"^ 




s I 


Saturation densitv of the leaf 


kg m 


ML 


E 


Transpiration rate 


kg m s-"' 




^1 


Convection coefficient 


9.14 J m"^s""^^^oc""^ 






Boundary layer resistance 
coemcienx 


1/2 -1 
s m 


1/2 -1 
T ' L 


L 


Latent heat of vaporization of 
water 


J kg"-^ 






Total radiation absorbed 




HL-^T-^ 


R 


Black body radiation 






r 

a 


Resistance to water diffusion 
across the boundary layer of air 


-1 ' 

s m 






Resistance to water diffusion 
within the leaf 


-1 

s m 




r.h. 


Relative humidity 




— 


T 


Temperature 




e 


T 

cl 


Air temperature 


oc 




Z 


Leaf temperature 




0 


V 


Wind speed 


-1 

m s 




w 


Leaf dimension perpendicular to 
the wind direction 


m 


L 


e 


Emissivity of the leaf surface 






0". 


Stefan-Boltzmann constant 


5.67x10"^ W m"^ 


-2 -1 -4 
HL T -^0 ^ 



L = length 

M = mass 

T = time 

0 = temperature 

H = ML^T~^ 
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USER'S GUIDE FOR PROGRAM TRANS 



Identification 

TRANS - A Program which Displays the Relationship between Leaf 

Temperature and Transpiration Rate 

Authors - Suzanne Miller, Larry Gales, Peter Farnum, Graham Carey, 
Center for Quantitative Science in Forestry, Fisheries 
and Wildlife, University of Washington, Seattle. March, 1977. 

Purpose 

Program TRANS is the computer supplement to the module "Trans- 
piration and Leaf Temperature" by David M. Gates (which is based on Gates, 
1968). This paper stresses the relationship, through the leaf energy budget, 
between the transpiration rate and temperature of a leaf and the independent 
variables of the climate and plant (radiation absorbed, air temperature, 
internal diffusive resistance, leaf dimension, relative hiimidity, and 
wind speed) • The computer program is based on this analysis and enables 
the student to create figures similar to those in the paper for any combi- 
nation of independent variables he wishes to specify . Thus the paper 
gives the student a theoretical understanding of the relationship between 
leaf temperature and transpiration and the associavad program helps him 
assimilate that knowledge by applying it to particu'^ar situations of 
interest • 

Operation 

The user specifies values fo: five -v r: seven independent 
variables and ranges for the other two (called changeable variables). 
The value of the first changeable variable is fixed at the lower end 
of the specified range while the value of the other changeable variable 
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changes continuously over its range. For each point in the independent 
variable space the energy budget equation (//7 in the module) is solved 
for leaf temperature using Newton's method, and then the transpiration 
rate is computed using this leaf temperature (from the last term in 
Equation (7)). This process generates a curve of transpiration versus leaf 
temperature for the five constant independent variables and the minimum 
of the first changeable variable as the value of the second changeable 
variable changes over its range. The value of the first changeable 
variable is then incremented to its minimum plus 1/3 of the range and 
the process is repeated to generate a second curve. 

This process is repeated to generate two families of four 
curves each of leaf temperature versus transpiration rate. For both 
families one changeable variable is fixed at its minimum, its minimum 
plus 1/3 of its range, its minimum plus 2/3 of its range, and its 
maximum value while the value of the other changeable variable changes 
continuously over its range. 

Example: let the two changeable independent variables 
be air temperature and wind speed and let their specified 
ranges be 10 to 40**C and 0.1 to 3.1 meters sec"^ . Then 
one family of four curves consists of transpiration versus 
leaf temperature for air temperature equal to 10, 20, 30 
and 40 degrees while wind speed varies continuously from 
0.1 to 3.1 meters sec"*^ . The other family contains four 
curvec for wind spited equal to 0.1, 1.1, 2.1, 3.1 meters 
sec"*^ while air tempvirature varies continuously from 10 
to 40**C. 
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Program Organization 

The program is organized according to the following flow chart: 



- MNODFLT " .T,?) 
I No 



Yes 



Read default values for 
all variables from the 
built-in default file 



Read in the next 
user-supplied data set 



( FINIS « > 
I No 



Yes 



Terminate 
Program 



Output 
appropriate 
error 
messages 



Check for errors in 
the input set just read 



i 

Yes ' w 
( Were errors found? ) 

\ No 



Calculate all the points 
in the leaf transpiration 
diagram. The points should 

all fall into 8 lines: 
4 lines where VAR(l) takes 
on 4 equi-spaced values, 
and 4 lines where VAR(2) 
takes on 4 equi-spaced 
values. Also, flag any 
failures to converge 



Write out the x,y,z 
coordinates for the 8 lines 
(z takes on the values 1-8) 
plus the titles for the plot 



Call the printer plot 
subroutine QQPR3D 
which generate^ the Plot 
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Input 

All input is handled by the format free input package (Gales and 
Anderson, 1978) which permits a user to assign values to variables by a 
."name^value" convention. Not all variables need be explicitly assigned by 
the user, however, as unassigned variables automatically assume default 
values. The input consists of any number of data sets, each of which is 
terminated by a dollar sign ($), Each data set generates a separate printer 
plot. 

The input for TRANS is divided into three classes: (a) variables 
having biological significance: VAR, VARMIN, VARMAX, QA, TA, RH, RL, D, W, 
and V; (b) variables which control certain program operations, such as 
program termination or the handling of default input: IPRINT, ECHO, NODFLT, 
and FINIS; and (c) variables which control the printer plots (default values 
are in parentheses): XMIN (0), XMAX^O) , YMIN (0), YMAX (0), ZMIN (0), 
ZMAX (9), XRICH (0.1), YRICH (8E-6), DFAULT (0), OVPRNT (.F.), AVE (.F.), 
INT2D (.F.), NX (60), NY (A5), and ZMAP (0,1,2,3,4,5,6,7,8,9). The variables 
in the first two classes are explained in the following INPUT TABLE, whereas 
the printer plot variables are explained in the user's guide for PRNT3D 
(Gales, 1978). 
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INPUT TABLE 



Name 


Type and 
Dimensions 


- Range Limits 


Description 


- VAR 


Integer (2) 


1,7 


Number codes for the first 
and second changeable 
variables. Default vaIupr 
are: 

VAR(l)-4, VAR(2)»7 


VARMIN 


Real (2) 


Depends on VAR 


Minimum values for the 
first and second change- 
able variables specified 
by VAR (1) and VAR(2). 
Default values are: 
VARMIN(l)-0, VARMIN(2)«=0.1 


VARMAX 


Real (2) 


Depends on VAR 


Maximum values for the 
first and second change- 
able variables specified 
bv VARQ) and VARf2^ 
Default values are: 
VARMAX(1)»1002, 
VARMAX(2)-5 


QA 


Real 


200,1200 


RflHlpnt! ATIATOV flViem*Vior1 

by leaf, in watts/meter^. 
Default value Is: QA-815. 
(Note: the code number 
for QA Is 1.) 


TA 


Real 


0,50 


Air temperature (^C). 
Default value Is: TA«40 
(Note: the code number 
for TA Is 2.) 


RH 

(-#3) 


Real 


0.1 


Relative humidity (as a 
fraction) . Default value 
Is: RH-0.5. (Note: the 
cuae numuer zor Kn XS j» ) 


RL 

(-#4) 


Real 


0,10000 


Diffusive leaf resistance, 
In sec meter" ^. Default 
value not specified as RL 
Is a changeable variable 
by default. (Note: the 
code number for RL Is 4.) 
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Input Table (cont.) 



Name 


Type and 
Dimensions 


Range Limits 


Description 


D 


Real 


0 < D <_ 2 


Leaf dimension parallel 
to wind in meters. 
Deiaux u value is : u*u • kjd 
(Note: the code number 
for D is 5.) 


W 

(-#6) 


Real 


0,10 


Leaf dimension transverse 
to wind, in meters. 
Default value is: W-0.05 
VNote: the code nimiDer 
for W is 6.) 


V 

(-■'■) 


Real 


0 < V _< 10 


Wind^ speed, in meters 
sec" . Default value is 
not specified as V is a 
changeable variable by 
default. (Note: the 
cone numoer ror v is / • j 


IPRINT 


Logical 


• F « y • T • 


A Logical value which 
causes the current values 
for all input variables 
(default as well as 
current user input) to 
oe prinLea* Deiauit value 
is: IPRINT « .F. 


ECHO 


Logical 


• F • } • 7 • 


A Logical value which 
causes the user's input 
to be echoed if ECHO«.T., 
or suppresses echoing if 
ECHO«.F. Default value 
is: ECHO-.T. 


NODFLT 


Logical 


• F « J • T • 


A Logical value which 
suppresses the input of 
default values if 
NODFLT-.T. Default value 
is: NODFLT-.F. 


FINIS 


Lonlcal 


• F • y • X • 


A Logical value which 
causes program termination 
if and only if FINIS«.T. 
Default value is: 
FINIS-. F. 
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The last four variables deserve special mention. 
The logical variable IPRINT controls the output of all Input variables 
which are currently In effect (default values as veil as those specified 
In the current Input set) • Setting IPRINT-TRUE (or T or .T.) displays 
the Input variables; setting IPRINT-FALSE (or F or .F.) suppresses the 
display. 

The logical variable ECHO controls the echoing of the Input cards. 
Setting ECHO=TRUE causes the subsequent Input set to the echoed; 
setting ECHO«FALSE suppresses the echo for the subsequent Input set. 
The logical variable NODFLT can be used to inhibit the automatic 
assignment of default values to input variables. If NODFLT is set 
TRUE in the current input set, then the current input set is assigned 
default values as usual, but all subsequent input sets merely 
accumulate more input values. In effect, the input values which exist 
after the 1-th input set is read, become the default values for the 
(l+l)~th input set. The standard default values may then be restored 
by setting NODFLT«FALSE, but, again, the effects of this change are 
d#;.layed until the next input set is read. To a limited extent, NODFLT 
permits a user to set up his own default values and can be very useful 
for executing a number of input sets which differ only in a few para- 
meters. Consider the following example in which a user wishes to view 
the same plot in a variety of ways by "zooming in" on various parts 
of it. Since the calculations are the same for all plots, only the 
viewing coordinates need be changed; 
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/INPUT SET 1: THE FOLLOWING VALUES BECOME THE DE FACTO/ 
/DEFAULTS FOR ALL SUBSEQUENT INPUT SETS: / 
NODFLT - TRUE, XRICH - 0.05, YRICH - 2E-6, 

VAR - 4, 7, QA - 700, TA - 30, RH - 0.2, D - 0.05. W - 0 05 
VARMIN - 0, 0.1, VARMAX - 6.1, 600, $ ' 

/INPUT SET 2: ZOOM IN (SEE RUN NO. 2 IN SAMPLE RUNS):/ 

XMIN - 27.3, XMAX - 33.3, YMIN - 4.296E-5, YMAX - 21E-5, $ 
/INPUT SET 3: ZOOM IN ON ANOTHER PART OF THE PLOT:/ 

XMIN - 20, XMAX - 27, YMIN - 8E-5, YMAX - 14E-5, $ 
/INPUT SET 4: ZOOM IN ON YET ANOTHER PART:/ 

XMIN - 25, XMAX - 28, YMIN - 40E-5, YMAX - 50E-5, $ 
/INPUT SET 5: STOP/ 

FINIS - TRUE, $ 

4, The logical variable FINIS controls program termination. The user 
should add the card: 

FINIS « TRUE, $ 

as the very last Input set. If FINIS is not set, the program will 
terminate abnormally. 

Output 

TRANS produces sets of plots, via subroutine PRNT3D, which are 
similar to those in the paper by Gates. Each plot contains a title, legend. 
X and y axis annotation, and printer plot lines. The title displays the 
values, and ranges of values, for the variables used to generate a plot and 
is best explained by an example. Consider the plot in the first sample 
run, whose title reads: 

QA-700.000 TA - 30.000 RH-. 200 D-.050 W-0.50 
***** THE ABOVE PARAMETERS ARE HELD CONSTANT***** 
LINES 1-4: RL-.OO 200.00 400.00 600.00 V».10 TO 6 10 
LINES 5-8: V-.IO 2,10 4,10 6.10 RL-.OO TO 600.00 ' 

ERIC 
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The plot contains eight lines with line 1 indicated by a string of I's, 
line 2 by a string of ... , and line 8 by a string of 8's. 

For this particular plot, the variables QA, TA, RH, D, and W are 
held constant at their stated values for all eight lines. The values for 
RL and V, however, are allowed to vary. For lines 1 through 4, RL is 
held fixed at a given value (0 for line 1, 200 for line 2, 400 for line 3, 
and 600 for line 4) while V varies continuously from 0.1 to 6.1. For 
example, line 3 features RL fixed at 400 while V varies from 0.1 to 6.1. 
For lines 5 through 8, the situation is reversed, and V is held constant 
at a given value (0.1 for line 5, 2.1 for line 6, 4.1 for line 7, and 
6.i for line 8) while RL varies continuously from 0 to 600. The user can 
select which variables are held fixed and which are allowed to vary. 

The plot legend, in conjunction with the numbers along the x and y 
axes, allows users to interpret the plot numerically. For example, the 
first line of the plot legend for the first sample run (page 51) reads: 

SCALE FACTORS « X-AXIS :E+01, Y-AXIS :E-04, Z-AXIS:E+00 
hence the point (x»2.133, y-4.707) [this is among a group of points marked 
with an "8" near the upper left corner of the plot] is interpreted as 
(x-2.133 X 10^ y»4.707 x 10"^) » (21.33"C, 0.0004907 kg/m^-sec). The next 
two lines in the plot legend indicate the number ("-9" means "more than 99") 
of points mapped to each z-axis level. 

The plot in run 1 also illustrates the problem of low resolution » 
which often plagues printer plots. The lines near the lower tight corner 
of the plot are so crowded together that it is difficult to distinguish 
- them (for example, line 2 contains only one, ^/isible point) . The "zoom in" 
feature of PRNT3D is invaluable here in that it permits users to blow up 
selected regions of the plot into as much detail as is desired. For example, 
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In run 2 the user input the coordinates 

(XMIN-27.3, YMIN-4. 296x10-') . (XMAX-33.3, YMAX-2l3clO-') 
of a rectangular window which enclosed the region he wished to see, and 
the computer responded with a full page blow up of the specified region. 
Zpom in and other features of PRNT3D are discussed in more detail in its 
user's guide. 

Restrictions 

The user musi restrict all inpuc variables according to the range 
limits l:isted in the input table, in o-ler to avoid unrealistic physical 
values. 

Error Message s 

There are four types of errors which may occur when attempting 
to execute program TRANS: 

1. Syntax errors in the user's input 

2. Range check errors 

3. Convergence errors 

4. Plot parameter errors 

For type 1 and 2 errors, the program flags the error, skips the calculations 
and plotting, and then reads the next input set. For type 3 errors, the 
program continues with the calculation and plotting after outputting an 
error message. For type 4 errors, the program suppresses plotting, outputs 
the error message, and reads the next data set. For a complete description 
of type 1 and type 4 error messages and actions, refer tc ::he user's guides 
for the format free input package and printer plot package, respectively. 

Range check errors occur If the range limits set for input 
variables are exceeded. These messages are of the form: 
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-™ ERROR NO, X IN SUB, INCHK 

yy,..y OUT OF RANGE 
yy. . .y ■ dd. . .d 

where x is an error number ranging from 1 to 11, yy...y describes the 
\iariable in question, and dd...d is the input value. The correspondence 
between error numbers and variables is as follows: 



Error No. 


Associated Variable 


1 


QA 


2 


TA 


3 


RH 


4 


RL 


5 


D 


6 


W 


7 


V 


8 


VAR(l) or VAR(2) 


9 


VARMIN(l) or VARMIN(2) 


10 


VARMAX(l) or VARMAX(2) 


11 


VAR(l) and VAR(2) 



Convergence errors occur when the values of combinations of input 
variables, which escaped the range checks, cause the iterative solution for 
leaf temperature to fail to converge. These messages are of the form: 

ERROR NO. 1 IN SUB. TLEAF 

CONVERGENCE REQUIRES MORE THAN d...d ITERATIONS 
QA » dd. . .d 
V » dd.,.d 
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LEAF TEMPERATURE - dd, ..,d 
NEWTON-S INCREMENT - dd...d 

ERROR NO. 2 IN SUB. TLEAF 

SLOPE - ZERO AT ITERATE - dd...d 

Sample Runs 

The annotated listing starting on the next page illustrates the 
control cards and input cards for two sample runs. Each input set is 
terminated by a $ and generates a one-page printer plot. These plots 
appear on the next few pages. 
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XTRANS»CM70000»T20. LEAF TRANSPIRATION 

ACf OUNT»XXXXXXXX,XXXXXXX,B, 

C0»1MENT. 

COMMENT. ***************************************************** *********** 

COMMENT,* THE FIRST CARD ABQVF IDENTIFIES THE JOB» SPFCIFIES THF ♦ 
COMMENT.* MEMORY PEOU IREMtNTS (55000 0CTAL)» AND ESTIMATED CENTRAL ♦ 
CO'^MENT.* PROCFSSQR TIMt IN SECONDi (20). THE SECOND CARD lOPNTIFlES * 
CO'^MEMT.* THE ACCOUNT AND PASSWORD. ♦ 
COMMENT .************* *************************************************** 
COMMENT. 

ATTACH»BPR3D»ID»BPR3D . 
ATTACH»BPF»ID»BFF. 
ATTACH»BTRANS. 
COMMENT. 

^2n**,^HJ ************** *************************************************** 

COM-iENT.* THF ABOVE CARDS ATTACH THE PRINTER PLOT ROUTINE. BPR3D» THF ♦ 
COMMENT.* CRCt PiDftM IMPUT ROUTINE* BFF» AND THE LEAF TP AN?! PI o ATI ON ♦ 
COMMENT.* ROUTINE* BTRANS. THEY APE ALL IN BINARY ^ORft * 
CO"^ ME NT. ******************************************** ******************** 

COMMENT. 

L0AD»«TRANS»3PR3n,BFF . 

EXECUTE. 

C HMMENT. 

comment!* THE ABOVE CARDS LOAD THE ROUTINES INTO MEMOPY AND CAUrE ♦ 
COMMENT.* CO^'TRCL TO Be PASSED TO TPANS POR EXECUTION. * 
CUMME NT. ******** ******************************************************** 
CO^^MENT. 
*EnR 
/ 

THE FOLLOWING TWO PLOTS A^E SIMILAR T'l THE PLOTS 
GATES (1966), ALTHJUGH THE UNITS USED THE»E WERE 
HERE, WE USE MKS. THE SEC^^O PLOT IS A DETAILED 
FIRST PLOT. 



FOUND IN 
CGS WHFREAS 
VIEW n«= THE 



THE FOLLCWING DEFAULT VALUES ARb ASSUMED FOP EAPH 0^ THE PLOTS 
UNLEiS rVERPIDDFN Y INPUT VARIABLES 

VAP(l) » V4R(2) » 7, 

QA » P15.8f, TA » tiOf RH « 0.5, M « 0, D » 0.05, W « 0.05, 
V » .000001, IPRINT = .F., 
VAP^TN » 0,0.1, VARMAX « 1002,5, 
ECHO » .T., FINIS » .F., KHDFLT 



0, ZM»y 
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/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 

/«4'******4***********************RUN 
/ 

/ THF FOLLOWING VALUES OVERRIDE THE DFFAULTS FOR THF FIRST RI'N 
/ 

OA » 700, TA » 30, RH « 0.2, VARMAX « 600,6.1, 
XRICH » 0.5, YRICH » 2E-6, « 

/ 

/«««4r«««*««««««««««*«i»«*«<f**«««**RUN 
/ 

/ THE FOLLOWING VALUES OVERRIDE THE DEFAULTS FOR THE SFCONO RUN 

6n 



NX » 60, NY » *t5f ZMAO . C, 1 , 2, 3, i , 6 , 7, 8 , Q 
XMIN « 0, XMAX « 0, YMIN = 0, YMAX « 0, ZMIN 
XRICH « 0.1» YRICH « 15-6, DFAULT » 0, 
nVPRNT » .F., AVE « .F., INT2D » .F. 



/ 
/ 

/ 

/ 
/ 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
*l 
I 
I 
I 



f 

*/ 
/ 
/ 

/ 
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lORINT • .T., 

OA ■ 700» TA ■ 3J» RH ■ 0,2f 

VARHAX - 600*6. If 

XBICH ■ 0.03» YRICH ■ lE-6t 

XHIN ■ 27. ^# XMAX ■ 33.3* YMIN ■ 4.296E-5# YMAX ■ 21E-5t 4 




FINIS ■ .T.# S 

♦ EOR 

♦ EOF 
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PRDe«AH -TRANS- READY FOR INPUT 

THE FOILOHING TWO PLOTS AR? SIMILAR TO THF PLOTS FQUNO, IN ^ 
GATES (l<)68>» ALTHOUGH THE UNITS USED THFRE WERE CGS MHFPEAS 
MERE* WE USE MK5. THE SECOND PLOT IS A DETAILED VIEW THE 
FIRST PLOT. 

THE FOLLOWING DEFAULT VALUED ARF ASSUMED FOR EACH Qp THE PLOTS 
UNLESS OVERRIDDEN BY IN«>UT VARIABLES 



VARd) ■ A, VAR(2) ■ 7» 

OA ■ eib.«e» TA ■ 40» RH • 0.5» RL ■ 0» D » ( 

V ■ .000001» TPRTNT ■ .F., 

VARMIN ■ 0»0.1» VARMAX ■ 1002»5» 

ECHO » .T,, PINTS ■ .F.» NOOFLT ■ .F.» 

NIX ■ 60» NY ■ A5» ZMAP ■ 0.» 1» 2» 3» <»» 5» fc »7» 8»9 

XMIN ■ 0» XMAX « 0» YMIN » 0» YMAX ■ 0» ZMIN 

XRICH ■ 0.1» YRICH ■ lE-6» OFAULT ■ 0» 

HVPRNT ■ .F., AVE » .F.» INT2D ■ .F. 



m 0, 7HAX ■ 0> 



^4 ********* ***9[iH I********** ******^> ********* ******* 

THE FOLLOWING VALUES OVFPRIPE THE DtFAULTS FOR THE CIPST 9UN 



Q4 • 700f TA ■ 30» RH ■ 0.2» VARMAX 
XRICH ■ 0.5* YRICH » 2E-6» S 



■ b00»(>,l9 
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QA" 700.000# TA> 30.000# RH« •ZOOr 0 • .050# W > ,0^0, 

THE ABOVE PARAMETERS ARE HELD CONSTANT ♦♦♦♦♦ 

hiNES l-'iiRL- 0 200.00 .<»00.00 600.00 V ■ .10 TO C.IO 

LINES ■ .10 2.10 • ^.10 6.10 RL» 0 TO bOO.OO 



L 
E 
A 
F 

T 
R 
A 

N 
S 
P 
I 
R 
A 
T 
I 
0 
N 



• 

4> 



1.662 
X— 
,966Y8 
18 
1 
1 



2.133 
• — X— • 



2.<i3<i 
• — X — 



2.735 
• — X — 



3.036 
X 



3.337 
X 



3.638 
X 



5.336 



^.707 



^1.078 



( 3.<i<i9 



K 
G 
/ 

< 2. 



620 



2. 191 



S 
E 
C 
G 

N 1. 

0 
) 
) 
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933 



,<i30Y 

X — 
1.862 



8 

88 



1 
1 
1 

77 
177 



8 

6t3 
SU 
88 
88 
88 



I 77 36 

II 777 88 

1 77 08 
1 77 8fc 
1 77 88 

I 77 86 

II 77 83 
6 7 88 



666 

llbo 
1 66 
1 6b 
11 
1 

11 
1 

11 
1 



o6 
666 
66 



77 88 
77 83 
77 88 
77 88 
77 OB 
7 88 
77 8b 
666 77 88 
66 77 8 
66 77 8 
1 660 77 8 

I 66 77 8 

II 666 788 
Hi 66 788 

111 66 788 

111 66 loo 

111 06788 
ill 66788 
ill 6786 
&i)5i>5 678d 
55555 668 
555580 

885555 
686 35555 
886333333555555 
8866 555555 
87o6<t'f'i ^'i ^ 

X X " — X— — — X— X — — X 

2. 133 2.'i3'i 2.735 3.036 3. 337 3. 636 

LEAr TEMPERATURE {uEGREES CcNTiGRADE) 



SCALE FACTORS ■ X-AXIS: E + 01 Y-AXISJ E-O^t Z-AXxi: c + oO 
Z0-2^« Ot-9)i> 1.000( 52)# 2.000( 0)> 3.000( o)» t.000{ 6) 

Z5-Z9- 5.000(35)/ 6.000(51), 7.000(52)/ 8.000(81)/ 9.000( 0) 
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PROGRAM -TRANS- KEAOY fOR INPUT 

4>»«* «*RUN Z* **************** *************** 

/ THE FOLLOWING VALUES OVERRIDE THE DEFAULTS FOR THE SECOND RUN - , 
/ 

IPRINT - .T.> 

OA ■ 700# TA " 30» RH ■ 0.2* 

VARhAX - 600#6.1# 

XRICH ■ 0.03* YRICH ■ lE-6* 

XHIN ■ <J7.3# XMAX ■ 33. 3# YMIN ■ <r.£96E-5# YMiX ■. 21E-5* X 
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THE COMPLETE LIST OF THE CURRENT VALUES OF INPUT VARUdLES FDR TRANS 



OA 

TA 

RH 

RL 

0 

W 

V 

VARHIN 

VARHAX 

ECHO 

FINIS 

NGOFLT 

NX 

NY 

ZMAP 



XMIN 

XMAX 

YMI N 

YMAX 

2MIN 

2KAX 

XRICH 

YRICH 

DFAULT 

OVPkNT 

AVE 

INT2D 

IPRINT 

% 



.7000000£*03# 
.3OC00OOL*0^# 
.2000000E*00# 

.^000000t-0i# 
•3C0O0O0E-01# 
• i00uCl)0E-03# 
0 / 
.6000000E^03# 

T* 

ft 

F# 

60# 
6, 

.27300001+02/ 
.3330OO0t^02/ 
.<»296000E-0<., 
•2i00000k-03# 

.900J000E+01/ 
• 3000000fc-01/ 
•lOOOOOOE-O^/ 
0 / 

F# 
F# 
F# 
It 



. 1000000E-*-00# 
.6100000E*01# 



It 
9t 



It 
6# 



3# 
It 
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Qkm 7O0.000# TA" 30.000# RH" .200* D ■ .0?0# W • .050* 

THE ABCiVE PARAMETERS ARE HELD CONSTANT ,^ ^ 

IHb Abuvt r* 200.00 400.00 600.00 V • .10 TO 6.10 

.10 2.10 " 4.10 6.10 kL" 0 TO 600.00 



LINES 1-4IRL- 
LINES 5-8IV 



L 
E 
A 
F 

T 

R 

A 

N 

5 

P 

I 

K 

A 

T 

I 

0 

N 



K 
6 
/ 
t 

.1 
♦ 

S 
E 
C 
0 
N 
0 
J 
) 



2.730 2.822 2^.923 3.025 3.127 

X X X X X— 

2.100Y 88 
I 8 
I 8 
177 68 
I 77 08 
77 88 
77 88 
77 8 
77 8 
77 8 
77 do 
77 88 
6d 77 80 
66 77 88 
6 77 8 
Db 77 8 
o6 77 6 
66 7768 
66 7788 
d66 7788 
66 778 
66 778 
D678fa 
6676tl 
667'J8 



i.228 
• — X— 



3.330 
— X 



1.910Y 
1 

' I 
I 

1.720Y66 
I 66 
I 
I 
I 

1.531Y 
1 



1 

1.341Y 
I 

lb 

I 51>t>55 



b3t;55 6788 
5b55b5 6788 
^bbb55 fc76 
5555tj78 

58855 

8&65355b 
686 2^255tj5> 
6666 2555t>3i» 
8766 2535355 

87d6 ^^^^^^^ct- 
837 6 55555 . 

8o7 66 '5 
087 66 
887 66 
8 /7336633 J3i3 33333333333 3 

877 66 
8877 66 
88 7 66 
88 7 6o «,«,c,«,<r<r44t't 

.OOY ^ ^lllt'ltltlit X 

7 llo z'lz'l 2T923 3T025 " 3.127 3.228 3.330 

LEAF TtMPERATORt (DEGREbS CtNlKiKADE) . 



1.151Y 
1 
I 
I 
I 

.961Y 
I 
I 
I 
I 

. 771Y 

I 
I 
I 

.581Y 
I 
I 
I 



SCALE FACTORS ■ X-AXiS: t+01 
^0-Z4- 0(-9)/ l.000( 0), 

Z»/-Z9" 5.000 (d5)/ 6.000(59)/ 



Y-AXIS« £-04 Z-AXISJ E+UO 
2.000( 5)/ 3.000120 ), 'i. 000( 18 
7.000(61)/ 8. 000173)* 9.000( 0) 




FINIS • 
06RA1 -T9ANS- TERMINATED 
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